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Subject:  Effect  of  Turbulent  Boundary  Layer  Flow  on  Measurement  of 

Acoustic  Pressure  and  Intensity  (Revised) 


Preface:  Because  of  an  error,  ARL/PSU  TM  83-163  has  been  revised. 

This  error  is  significant  enough  to  alter  the  conclusions,  so  it  is 
necessary  to  re-issue  a  revised  version  of  the  report  which  is  contained 
herein.  Please  discard  ARL/PSU  TM  83-163,  and  use  this  version  in  its 
place. 


Abstract:  Acoustic  measurements  performed  with  pressure  sensors 

(microphones  and  hydrophones)  can  be  subject  to  misinterpretation  if  the 
sensors  are  in  contact  with  a  turbulent  boundary  layer  flow.  This 
misinterpretation  is  a  result  of  non-propagating  pressure  fields  (flow 
noise)  generated  by  the  turbulent  motions  of  the  fluid.  In  this  paper, 
a  simple  analysis  is  given  whereby  the  flow  noise  response  of  a  pressure 
sensor  placed  in  a  turbulent  boundary  layer  flow  can  be  estimated.  If 
the  purpose  of  the  sensor  is  to  measure  the  sound  emitted  from  an 
independent  source  outside  of  the  turbulent  boundary  layer,  then  a  bias 
for  the  measurement  can  be  calculated.  A  bias  formula  for  intensity 
measurements  made  with  two-sensors  under  a  turbulent  boundary  layer  is 
also  derived  under  the  assumption  of  a  very  low  Mach  number  flow. 
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INTRODUCTION 

There  has  been  considerable  interest  recently  in  the  measurement  of 
acoustic  pressure  and  acoustic  intensity  in  ducting  systems  and  other  systems 
that  sustain  a  mean  flow  [1-5]*.  The  work  of  Munro  and  Ingard  [3]  is 
fundamental  and  provides  a  mathematical  description  of  acoustic  intensity  for 
a  medium  in  motion.  They  point  out  that  two  potential  problems  exist  when  one 
attempts  to  measure  intensity  in  a  system  with  flowing  fluid.  The  first  deals 
with  the  definition  of  intensity,  where  the  majority  of  their  paper  addresses 
this  issue.  The  second  problem  deals  with  what,  exactly,  does  a  pressure 
sensor  (microphone  or  hydrophone)  which  is  placed  in  a  flow  field  measure? 

This  problem  is  associated  with  "wind  noise"  or  "flow  noise"  and  is  not 
discussed  further  in  Ref.  [3]. 

For  the  experimentalist,  flow  noise  can  be  an  important  problem  that 
deserves  some  attention.  It  must  be  pointed  out  that  it  is  a  difficult 
problem  to  solve,  in  general,  because  each  situation  is  likely  to  have  unique 
features  associated  with  the  flow  field,  and  with  the  acoustical  measurement 
apparatus.  In  an  attempt  to  generalize  the  latter,  we  can  usually  envision 
the  measurement  sensors  to  be  situated  either  in  the  free  stream,  or  somehow 
flush  mounted  in  the  walls  of  the  system.  For  example,  measurements  of  tire 
noise  may  be  conducted  on  a  moving  vehicle  with  a  microphone  placed  in  the 
free  stream  near  the  tire  and  roadway.  This  situation  is  discussed  by 
Oswald  and  Donavan  [4],  where  the  errors  due  to  the  probes  being  in  a  free 
stream  are  assessed.  They  conclude  that  very  good  accuracy  is  achieved  with 
intensity  probes  for  those  situations  in  which  the  flow  Mach  number  is  less 

*Numbers  in  brackets  denote  references  listed  at  the  end  of  the  paper. 
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than  0.1  and  the  turbulence  intensity  is  less  than  about  5  percent.  These 
results  are  based  on  anechoic  wind  tunnel  measurements  using  microphones 
equipped  with  streamlined  nose  cones. 

The  effect  of  turbulence-induced  noise  on  flush-mounted  sensors  has  been 
of  concern  to  those  involved  in  underwater  sonar  system  design  and  evaluation, 
and  to  researchers  that  utilize  wind  or  water  tunnel  test  facilities  in  aero 
and  hydroacoustic  investigations,  e.g.,  [6-8].  The  predominant  source  of  this 
noise  is  the  turbulent  boundary  layer  that  exists  between  the  walls  of  the 
system  and  the  free  stream.  Boundary  layer  turbulence  produces  pressure 
flucutations  that  span  wide  frequency  and  wavenumber  ranges.  The  response  of 
a  flush-mounted  sensor  placed  beneath  a  turbulent  boundary  layer  depends 
therefore  on  the  details  of  the  boundary  layer  wavenumber/frequencv  pressure 
spectrum  and  on  its  own  spatial  averaging  characteristics.  One  of  the  first 
mathematical  descriptions  of  this  type  of  interaction  is  given  by  Uberoi  and 
Kovasznay  [9]. 

Many  noise  control  engineering  problems  involve  the  measurement  of  an 
acoustical  phenomenon  associated  with  fluid  moving  devices  or  systems.  Such 
problems  include  the  diagnosis  and  reduction  of  noise  sources  of  aircraft, 
ships,  cooling/heating  systems,  and  hydraulic  systems,  to  name  a  few.  Other 
problems  are  associated  with  the  evaluation  of  materials  or  products  that 
are  to  be  subjected  to  flow.  This  latter  area  has  led  to  the  recent  interest 
in  duct-type  test  facilities  [1,5].  The  identification  of  boundary  layer- 
induced  transducer  response  would  be  of  interest  in  most  of  those  applications 
that  use  flush-mounted  sensors.  In  this  paper,  a  simple  analysis  of  this  type 


of  response  is  presented.  The  analysis  makes  use  of  research  findings  in 
several  disciplines:  hydrodynamics,  hydroacoustics,  and  acoustics.  Well- 
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known  concepts  of  each  discipline  are  combined  to  provide  a  useful  method  for 
estimating  a  bias  that  can  exist  in  such  measurements.  Numerical  examples  are 
presented  for  typical  transducer  sizes  and  air  flow  velocities  encountered  in 
pract ice . 

PRELIMINARIES 

The  problem  to  be  investigated  is  the  determination  of  that  portion  of  a 
pressure  or  intensity  spectrum  due  to  turbulent  boundary  layer  pressure 
fluctuations  when  the  measurement  is  performed  under  a  turbulent  boundary 
layer  with  flush  pressure  sensors  of  given  shape  and  size.  For  convenience, 
this  problem  is  approached  using  a  planar  surface  as  a  model.  The  results 
obtained  with  this  model  may  be  extended  to  duct  flow  problems  providing 
that  the  appropriate  changes  in  turbulent  boundary  layer  development  are 
accounted  for. 

Flow  Cons i der ations 

Consider  the  planar  flow  model  of  Figure  1.  The  mean  flow  velocity  is 
denoted  by  u^  and  it  is  assumed  that  u^/c  =  M  <  0.1,  where  c  is  the  sonic 
velocity.  The  fluid  (air,  water,  oil,  etc.)  is  characterized  by  its 
viscosity,  p,  and  density,  p.  The  kinematic  viscosity  is  v  =  p/p,  A 
parameter  which  helps  classify  the  dynamic  similarity  of  geometrically 
similar  flows  is  the  Reynolds  number,  Re^  =  Uy£/v.  Here,  £  is  a 
characteristic  length  of  the  plate,  tube,  or  body.  In  the  case  of  viscous 
flow  over  a  surface  (Figure  1),  a  shear  layer  develops  very  close  to  the 
surface.  In  traditional  descriptions  of  this  boundary  layer,  the  length  used 
in  the  Reynolds  number  may  be  either  the  streamwise  coordinate,  x^,  or  the 
thickness  of  the  layer,  6.  The  boundary  layer  is  defined  as  that  region 
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above  the  surface  where  the  velocity,  u  depends  on  x^ •  Its  thickness,  6,  is 

defined  at  that  value  of  where  u  =  0.99  Uq .  In  some  instances,  the  mean 

velocity  may  be  directed  at  some  angle  to  the  plane.  In  these  cases,  a  mean 

pressure  gradient  will  exist  along  the  surface  and  alternative  definitions  for 

the  boundary  layer  flow  are  necessary.  For  the  analysis  of  this  paper,  it  is 

assumed  that  the  flow  is  parallel  to  the  surface. 

Initially,  near  the  leading  edge,  the  boundary  layer  is  laminar.  The 

streamlines  are  straight  and  parallel  in  the  boundary  layer  and  no 

fluctuations  in  velocity  or  pressure  are  produced.  The  laminar  layer  persists 

for  only  short  distances  before  instabilities  arise  and  the  layer 

"transitions"  into  a  turbulent  layer.  The  transition  location  for  a  smooth 

flat  plate  can  be  estimated  from  the  transition  Reynolds  number,  typically 

5  6 

in  the  range  1 10]:  10  <  ReXt  <  10  .  In  the  case  of  a  smooth  pipe  flow, 

ReXt  =  2300,  where  the  velocity,  u^ ,  is  that  at  the  center  of  the  pipe.  The 

transition  location  is  sensitive  to  surface  finish.  In  general,  these 
transition  Reynolds  numbers  would  be  much  lower  for  rough  surfaces.  The 
turbulence  in  the  layer  is  assumed  homogeneous  and  of  zero  mean. 

Acoustic  Measurement  Considerations 

It  is  assumed  that  pressure  sensitive  transducers  are  mounted  flush  in  a 

surface  as  illustrated  in  Figure  1.  It  is  further  assumed  that  the  purpose  of 

these  sensors  is  to  measure  an  acoustic  field  (either  the  pressure  or  the 
intensity  vector),  and  that  this  field  is  statistically  independent  of  the 
pressure  field  generated  by  the  turbulent  layer.  The  rms  acoustic  pressure 
originating  from  the  source  of  interest  is  denoted  by  p  and  is  sensed  by  a 

d 

single  transducer.  The  turbulence  produces  rms  pressure  p^,  which  together 


..V 
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with  p 

a 


results  in  a  total  measured  pressure: 


P  =  Pa  +  PT  • 


i 

(1) 


The  acoustic  intensity,  I  ,  is  determined  in  the  usual  way  using  a  pair 

cl 

of  sensors  located  a  distance  Ax  apart  [11,12]: 

Ifl  =  ~  Im(Gl2}/pwAx  .  (2) 

Equation  (2)  contains  the  quad-spectrum  (imG.^)  as  measu:  .  jy  the  pair  of 
transducers,  and  to  is  the  radian  frequency.  It  is  assutr  '■t  all 

corrections  for  sensitivity  and  phase  mismatch  have  been  accounted  for. 


BIAS  DUE  _T0  TURBULENT  JOU N D A RY  LAYER  Jv LOW 

If  our  objective  is  to  measure  either  p  or  I  in  the  presence  of  a 

B  B 

turbulent  boundary  layer  (TBL) ,  Equation  (1)  shows  that  this  measurement 
may  be  biased.  If  6  is  a  measured  quantity,  and  0  is  the  desired  or  true 
quantity,  then  the  bias  is  defined  l 1 3 ]  as: 

b ( 9 ]  =  E{6}  -  9  ,  (3) 

where  E{  }  denotes  expectation  operation.  The  normalized  bias  is  b ( 9 j / 8 . 


Pressure  Spectr urn  _Me_a s urement 

The  one-sided  power  spectrum  measured  by  a  single  flush-mounted  sensor 
is  defined  [13]  as 


G 

PP 


(to)  =  lira  y  E{P* 


p} 


(4) 
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where  P  is  the  finite  Kourier  trnnsf jnn  of  p  and  the  asterick  denotes  complex 
conjugate’.  using  Equation  Cl)  in  Equation  (4)  results  in: 


g 

pp 


Gaa  +  GTT  + 


G  „  + 
a  r 


G_ 


I  a 


where  G ^  and  are  cross  spectra  between  the  TBL-induced  pressure  and  the 

acoustic  pressure.  Because  of  statistical  independence,  these  two  cross 
spectra  are  zero,  and 


G_  =  +  Gtt  • 

pp  aa  TT 


Co) 


An  estimate  for  the  acoustic  pressure  autospectrum  is  G  (w),  so  the  bias  in 
this  estimate  follows  directlv: 


b  [G  J  =  E{G  }  -  G 
'  aa ;  aa  aa 


=  G  +  G  _  -  G 
aa  TT  aa 


The  normalized  bias  is: 


‘  CTT  • 


VhJ  ‘  <WG 


(6) 


TT'  aa 


(7) 


Equation  (7)  is  an  expected  result  for  two  independent  sources;  it  states 
that  the  normalized  bias  error  for  the  acoustic  autospectrum  is  simply  a 
noise-to-signal  ratio,  where  the  spectrum  of  the  turbulent  boundary  layer 
pressure  field  represents  the  noise  term.  In  practice,  the  numerator  might 
possibly  be  determined  by  turning  off  the  acoustical  source  while  retaining 
the  flow.  When  this  is  not  possible,  G^^,  can  be  estimated  by  procedures 


i 

1 

i 

j 

\ 

! 
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given  in  Appendix  A.  Likewise,  the  denominator  of  Equation  (7)  could  be 
determined  by  measuring  the  acoustic  autospectrum  in  the  absence  of  flow. 
When  this  is  not  possible,  the  normalized  bias  for  G  can  still  be  computed 

cici 

with  a  knowledge  of  G  and  Equation  (8),  i.e., 

PP 

£,(G  )  =  (C  /C,  )/[l  -  G  /C  I  .  (8) 

b  *■  aa '  '  TT  pp  IT  PP 


Intensity  Spectrum  Measurement 

Letting  I  ( 6j)  be  the  intensity  spectrum  calculated  from  measurements 
made  with  the  pair  of  sensors  (labeled  1  and  2)  sketched  in  Figure  1,  it 
follows  that: 

b[I  1  =  -  ImfGpJ/pwAx  +  Iin(Gp  ]/' puiAx  .  (9) 

Here,  Gp  is  the  measured  cross  spectrum  between  the  two  sensors  with  flow 
and  Gp  would  be  the  cross  spectrum  measured  in  the  absence  of  flow.  Noting 
that 


P,  =  Pa  +  PT  UOa) 

1  al  ll 

and 

P,  =  Pa  +  PT  ,  OOb) 

a2  l2 

where  subscripts  1  and  2  refer  to  sensors  1  and  2  respectively,  we  find 
that 


12 


=  G 


12 


GI  t  + 
1  2 


V2 


+  G_ 


T1  a2 


i 


(11) 
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Again,  because  of  the  statistical  independence  between  p^  and  p^, 
Equation  (11)  reduces  to: 


12 


12 


+■  G, 


T.T 


1  2 


(12) 


If  the  two  sensors  are  placed  farther  apart  than  the  streamwise 
correlation  length.  A,  of  the  turbulent  pressure  fluctuations,  then  the 
turbulence-induced  cross  spectrum  G  is  nearlv  zero.  It  would  be 

Tr2 

concluded  then  [Equation  (9)]  that  the  estimate  for  the  acoustic  intensity 
is  unbiased  when  the  measurement  is  performed  under  a  TBE.  It  is  presumed 
that  the  measurements  of  Oswald  and  Donavan  [4]  were  performed  using 
probe  separations  greater  than  the  correlation  lengths  of  the  turbulent 
pressure  fluctuations  that  occur  in  the  free  stream  of  the  wind  tunnel 
used.  That  is  why  they  found  very  little  error  in  the  intensity 
measurements  even  when  the  noise-to-signal  ratio  for  the  noisiest  micro¬ 
phone  was  greater  than  unity.  Nonetheless,  it  is  important  to  examine 
Equation  (9)  in  some  detail  in  an  attempt  to  identify  those  ranges  of 
turbulence  length  scales  where  G  cannot  be  neglected. 

V2 

The  cross  spectrum  of  the  pressure  fluctuations  that  occur  at  the  wall 
under  a  TBL  was  modeled  by  Corcos  [6]  and  others.  For  this  discussion, 
the  Corcos  model  is  sufficiently  accurate  and  accounts  for  the  spatial 
averaging  effects  introduced  by  the  finite  size  of  the  measurement  sensors. 
This  model  states  that  the  measured  cross  spectrum  is  of  the  form: 


gt  t  5  GTT^)A^c^B,'kcn>e 


(13) 
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Here,  £  represents  the  separation  of  the  two  sensors  in  the  x  -direction 

and  n  is  the  separation  in  the  spanwise  [x^j  direction  (assumed  to  be  zero 

for  this  study).  The  wavenumber  k  is  u/u  and  is  called  the  convection 

c  c 

wavenumber  since  u  represents  the  speed  at  which  the  pressure  producing 
turbulence  is  "convected"  over  the  surface.  This  process  is,  in  general, 
dispersive  with  u^  being  a  weak  function  of  ui.  However,  the  dependence  is 
slight  and  many  workers  in  the  field  assume  (as  we  will  here)  that 
u£  =  0.7  u^ .  The  functions  A  and  8  are  normalized  and  can  be  approximated 
by: 


A(kc5)  =  exp(-  0.ll5|kcC| j 
B(kcn ;  =  exp(-  0 . 7  j k^n | ) 


(14) 

(15) 


With  the  intensity  probe  oriented 

thus,  the  imaginary  part  of  G„ 

T2 


in  the  direction 
is  of  the  form: 


of  flow. 


£  =  Ax  and  n  =  0 


Im(GTiT?)  =  ~  GXT(u)A(tccAx)sin 


(16) 


We  can  now  use  Equation  (12)  in  Equation  (9)  to  obtain: 


b(l  )  =  -  lm(G  ) / pwAx 
a  V2 

-  [GTT(u}A(kcAx)sin  k£  Ax  ]  /  pcoAx 


or 
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G  ( to)  sin  k  Ax 

-^T-~  A(k  <L.. 

pu  c  k.  Ax 

c  c 


(17) 


A  non-dimensional ized  plot  of  this  bias  is  given  in  Figure  2  where  k  Ax  is  the 

c 

independent  variabLe.  A  second  axis.  Ax/ A  is  also  shown,  where  the 
correlation  length  is  calculated  from: 


OO 

A  =  J  Afk  Fjdc 
0 


A  =  1 1 n/4k  ,  .  ( 18) 


(The  accepted  form  of  the  correlation  length  is  based  on  the  magnitude  of  the 
cross  spectrum.)  This  plot  shows  clearly  the  rapid  decorrelation  of  turbulent 
energy  as  Ax/ A  increases.  From  it  we  can  conclude  that  an  intensity  spectrum 
estimate  would  be  essentially  unbiased  for  sensor  separations  greater  than 
about  one  TBL  correlation  length  (A).  Because  of  Equation  (18),  we  must 
recognize  that  A  is  inversely  related  to  frequency  such  that  a  bias  can  exist 
in  intensity  measurements  for  low  frequencies  (f^j,  i.e.. 


£l  <  |i  >- 


u  =0.7  u., j 

K  n  (I  ' 


Below  f^,  the  magnitude  of  this  bias  may  still  be  small  because  of  the 
multiplicative  factor  G  / pu^  of  Equation  (17).  This  magnitude  can  be 
qualified  conveniently  by  examining  the  normalized  bias. 
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For  the  assumed  propagating  acoustic  field, 

Xm (^12  ^  =  IG12 I Sln  kAx  ’ 


(20) 


where  k  is  the  sonic  wavenumber.  But,  the  finite  difference  criterion  for  the 
two-sensor  intensity  method  requires  kAx  «  1.  Therefore,  we  can  let 
sin  kAx  =  kAx.  Then,  the  normalized  bias  is: 


,h.)  * 


G  sin  k  Ax 

A^Ax)  k"*Ax 
’  1 L  1  c 


(21) 


where  is  introduced  as  the  convective  Mach  number  (u^/c).  From  Equations 
(17)  and  (21),  the  ordinate  of  Figure  2  is  also  interpreted  as: 


b(l  )pu  „  I G  | 

=  c  fr  )V,  Ly : 

CT.r  h^a-’c  G 

T  f  TT 


(22) 


As  an  example,  if  Ax/ A  =  1.0,  Equation  (22)  would  have  a  numerical  value 
of  order  0.01.  If  the  air  velocity  is  on  the  order  of  30  m/s,  where  M  =  0.1 
(above  this,  the  intensity  definition  used  becomes  invalid  [3])  and 


Mc  =  0.07,  we  see  that 


c  ( i  )  .  i  Jit. 


Clearly,  e  (i  )  is  only  14  percent  for  the  extreme  condition  of  G_ ,  =  | G  _  j . 
If  Ax  is  much  larger  than  A,  these  results  suggest  that  accurate  intensity 
measurements  can  be  performed  under  a  TBL  even  when  the  sensed  TBL  auto- 
spectral  density  is  higher  than  that  due  to  the  acoustic  source.  That  is, 
the  normalized  bias  for  G  [Equation  (7))  could  be  as  much  as  two  orders  of 
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magnitude  higher  than  the  normalized  bias  for  1  .  As  pointed  out  by 
Corcos  [6],  and  detailed  in  Appendix  A,  these  biases  can  be  reduced 
significantly  by  using  large  diameter  microphones  in  the  measurements. 

CONCLUSIONS 

In  this  paper,  an  analysis  is  given  for  a  potential  data  interpretation 
error  that  may  occur  in  acoustical  measurements  that  employ  pressure  sensors 
mounted  in  walls  of  systems  that  sustain  low-Mach  number  flow.  A  review  of 
recent  literature  reveals  that  measurements  of  this  type  are  conducted  in  the 
noise  control  engineering  field,  so  it  appears  that  an  analysis  of  this  type 
may  be  of  interest.  Under  the  assumptions  that  the  flow  Mach  number  is  low 
(so  that  the  zero-flow  definition  of  intensity  may  be  used),  and  that  the 
pressure  fluctuations  generated  by  the  turbulent  boundary  layer  flow  are 
statistically  independent  of  the  acoustical  pressure  under  investigation,  the 
normalized  bias  expected  in  the  measurement  of  the  acoustical  autospectrum  is 
shown  to  be  simply  a  noise-to-signal  ratio.  The  bias  in  the  intensity 
spectrum  is  shown  to  be  significantly  smaller  than  this  noise-to-signal  ratio 
provided  that  the  sensor  separation  is  greater  than  the  streamwise  correlation 
length  of  the  pressure  fluctuations  created  in  the  turbulent  boundary  layer. 

Methods  for  predicting  the  TBL  wall  pressure  spectrum  as  measured  by 
transducers  of  finite  size  are  given  in  an  Appendix  along  with  example 
calculations.  For  the  typical  air  flow  velocities  and  sensor  sizes 
encountered  in  practice,  these  example  calculations  show  that  the  measurement 
of  acoustical  pressures  that  are  lower  than  approximately  60  dB  re  20pPa/  Hz 
will  be  influenced  by  the  presence  of  a  TBL. 


1 


t 
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The  limitations  of  the  analysis,  in  addition  to  those  associated  with 
the  assumptions,  include  the  following: 

(1)  Vibration  of  the  surface  may  lead  to  sensor  response  not 
accounted  for  in  the  analysis.  This  would  be  most 
critical  for  liquid  flows  where  fluid  loading  is  non- 
negligible . 

(2)  Angle  of  attack  effects  may  require  modifications  in 
the  flow  equations  and  pressure  spectrum  analyses 

of  Appendix  A. 
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APPENDIX  A 

Estimation  of  G-pfCcj) 

Equations  (7),  (8),  (17)  and  (21)  all  contain  G^Xw)  which  is  the  auto¬ 
spectrum  of  the  TBL  pressure  fluctuations  as  measured  by  a  flush-mounted 
pressure  sensor  of  given  physical  size.  If  the  source  of  interest  can  be 
turned  off  while  still  retaining  the  TBL  flow,  then  G^(u)  can  be  measured 
directly;  otherwise,  it  must  be  estimated  from  theoretical  analysis.  This 
appendix  gives  such  an  analysis. 

The  wall  pressure  spectrum  is  calculated  from  methods  first  presented 
by  Uberoi  and  Kovasznay  [9].  In  particular: 

00  oo 

Gtt(u)  =  /  /  «(k1,k3,o)|H(kl,k3)|2dk1dk3  ,  (A-l) 

—  OO  — OO 

2 

where  |H(kj,k3)|  is  the  in-plane  wavenumber  response  function  of  the  sensor, 
and  $(kj,k3,w)  is  the  wavenumber/f requency  spectrum  of  the  TBL  wall  pressure 
fluctuations  at  the  given  Reynolds  number. 

The  sensor  spatial  response  function  depends  on  the  shape  of  the  sensor. 
Most  microphones  and  hydrophones  used  in  pressure  and  intensity  measurements 
are  circular.  For  circular  elements  of  active  diameter,  D  =  2R, 

| H( kR) | 2  =  [2J1(kR)/kR)  ]2  ,  (A-2) 

2  2  2 

where  k  *  k^  +  k3,  and  is  the  Bessel  function  of  order  one.  Figure  3 

2 

shows  the  variation  of  | H ( k  ^  R ) |  with  k^R.  Because  of  polar  symmetry,  this 
variation  is  identical  to  that  in  the  k3R  domain.  The  measured  response 
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function  shown  in  Figure  2  is  by  Farabee  and  Geib  [14].  If  rectangular 
elements  are  used  in  a  given  application  [6]  , 

|H(k.L.)|2  =  sin2(k.Li/2)/(k.L./2)2  ,  (A-3) 

i  -  1,3  , 

and  is  the  length  of  the  side  in  direction  x^. 

Over  the  past  twenty  years  or  more,  many  papers  have  been  published  on 
the  wall  pressure  spectrum  of  turbulent  boundary  layers.  Even  today, 
disagreement  exists  among  researchers  as  to  the  true  form  of  this  spectrum. 
Measurement  of  <t>(k^  ,k^ ,  u)  is  extremely  difficult  because  of  its  multi¬ 
dimensionality  and  its  wide  dynamic  range.  Likewise,  theoretical  modeling 
usually  requires  some  experimental  data  in  order  to  fix  proportionality 
constants  that  arise. 

Aside  from  these  difficulties,  the  hydroacoustics  community  generally 
agrees  that  'tfk^.k^.w)  is  a  symmetric,  gradually  decreasing  function  in  k^, 
and  varies  in  u  like  narrow-band  random  noise.  The  variation  of  ^(k^.k^.wj 
with  k^  (streamwise  direction)  is  more  complicated.  Figure  4  is  a  schematic 
description  of  4>(k^ ,0, w) .  A  sharp,  well-defined  peak  occurs  at  the  convection 
wavenumber.  The  peak  is  simply  explained  in  terms  of  pressure  pulses  created 
by  the  eddies  of  turbulent  flow  that  travel  at  speed  u^ .  Because  the  convective 
wavenumber  is  much  greater  than  the  sonic  wavenumber  (w/c)  by  a  factor  of 
0(Mc*),  energy  at  the  convective  peak  is  non-propagating.  The  only  energy 
that  propagates  is  that  for  which  k^  <  to/c.  The  energy  in  the  sonic  wave- 
number  portion  of  the  spectrum  is  very  small  in  comparison  to  that  in  the 
high,  convective  wavenumber  region. 
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A  so-called  low  wavenumber  region  exists  between  the  sonic  and  convection 
wavenumbers.  The  energy  in  this  region  is  considered  important  because  this 
band  of  wavenumbers  (w/c  <  <  m/u^'  usually  couples  well  with  flexural 

wavenumbers  of  real  structures.  (The  coupling  is  most  efficient  for  under¬ 
water  structures.)  Figure  4(a)  indicates  a  parabolic  curve  that  relates  the 
resonant  frequency  of  a  typical  flat  plate  with  wavenumber.  The  constant 
a  is  related  to  the  bending  stiffness  of  the  plate. 

Measurements  of  <j(kc  ,0 ,  ojj  have  been  quite  successful  in  recent  years. 
Inspection  of  Equations  ( A— 1 )  through  (A-3)  reveals  that  a  measured  spectrum 
would  be  dominated  by  the  energy  at  k^  =  kc  if  very  small  pressure  transducers 
were  used  to  make  the  measurement.  The  spectra  measured  by  Blake  [151  with 
pinhole  microphones  on  a  plane  surface  in  an  air  flow  are  quite  representative 
of  O(kc,0,uj). 

Measurements  in  ranges  where  k.  <  k  are  less  common  and  very  difficult 

1  c 

to  perform.  Wavenumber  filters  are  necessary  in  order  to  block  out  the 
contributions  from  the  sonic  and  convective  regions.  Limited  low  wavenumber 
pressure  spectra  are  given  by  Blake  and  Chase  [8]  and  by  Martin  and  Leehev 

[16]. 

In  an  attempt  to  formulate  a  model  of  $(kj ,k-j,aij ,  Maestrello  [17] 
measured  space-time  correlation  functions  of  TBL  pressure  fluctuations  with 
pairs  of  small  microphones  in  a  planar  surface.  These  functions  were  then 
Fourier  analyzed  in  both  space  and  time  to  obtain  the  wavenumber/ frequency 
spectrum.  The  Maestrello  model  is  characterized  by  a  low  wavenumber 
region  that  remains  flat  with  decreasing  k^.  In  a  more  recent  analysis. 
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Chase  (18]  argues  that  the  spectral  level  should  eventually  decrease  to  zero 

~> 

as  k^  ♦  0.  His  model  shows  a  k^-  dependence  tor  k^  <  and  is  given  by: 


$'vk,  ,k3 , 


-  •Mhyfc5  *  cdk 


^  .  1.2  ',,,-5 

1  +  k3^ 


(A-4) 


K.  =  [(a)  —  u  k  )~/h2v2]  +  (kf  +  k2  )  +  fb.fi)  2  ,  ( A—  5 ) 

1  L K  cl'  L*J  v  i  3J  viy 

i  =  M,T.  The  velocity,  ,  is  the  friction  velocity  and  is  related  to  the 

r  i  /  2 

wall  shear  stress  by  v^  =  't  /p)  .  For  a  smooth  flat  plate, 

t  =  0.029  pu2Ke“l/5  (A-6) 

w  0  Xj 

where  the  Reynolds  number  is  based  on  the  distance  back  from  the  beginning  of 
the  TBL .  The  TBL  thickness  at  x^  is  approximated  by: 

o  *  0.37  x,  Re"1/5  .  (A-7 ) 

1  Xj 

The  constants  of  Equations  (A-4)  and  (A-3)  need  to  be  selected  through 

comparison  of  predictions  with  appropriate  data.  Chase  considers  the 

Martin-Leehey  data  [16]  to  be  good  at  low  wavenumbers.  Also,  using  the 

Blake  data  for  k,  =  k  ,  the  following  constants  are  established: 

1  c 

»  hT  -  20.0  ,  CM  =  0.03  ,  CT  =  0.004  , 


bM  =  0.75  ,  and  bT  =  0.38 
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A  comparison  of  the  Chase  and  Maestrello  models  is  given  in  Figure  5. 
Clearly,  large  differences  are  apparent,  but  it  must  be  emphasized  that  both 
models  depend  to  fair  degree  on  the  accuracy  of  the  experimental  data  used 
in  their  development.  Because  the  Chase  model  makes  use  of  the  data  base 
obtained  at  very  low  Mach  numbers,  it  is  the  recommended  model  for  the 
analysis  of  this  paper. 

Equation  (A-l)  is  solved  by  numerical  integration  for  the  flow  conditions 

considered  by  Blake  [15]  in  his  measurements  of  $fk  ,0,u).  The  integrations 

k  c 

are  terminated  at  (-10  k  ,  10  k  1  for  both  variables.  The  microphone  diameter 
used  by  Blake  is  0.79  mm  (pinhole  cap  over  standard  B&K  4138  microphone),  and 
his  airflow  velocities  range  from  22.2  to  50  m/s.  A  comparison  of  Blake's 
data  for  22.5  m/s  (ReXl  =  4.33  *  1 0 ^ )  with  the  prediction  using  the  Chase 
model  for  <t>(kj  ,k^ ,  w)  is  given  in  Figure  6.  The  agreement  is  reasonably  good, 
except  for  low  frequencies,  where  differences  of  order  6  dB  occur. 

Calculations  at  other  velocities  compare  equally  well  with  the  experimental 
spectra. 

To  further  illustrate  the  types  of  sensor  response  expected  when  the 
sensor  is  placed  under  a  TBL,  Figure  7  is  shown.  Three  typical  Reynolds 
numbers  are  selected  which  would  correspond  to  air  velocities  (standard 
temperature  and  pressure)  of  6.10,  18.29,  and  24.38  m/s  at  a  streamwise 
sensor  location  of  x^  =  0.30  m.  Spectra,  for  four  microphones  of  diameter 
ranging  from  3.18  to  25.4  mm  are  shown.  These  predicted  spectra  are 
essentially  flat  up  to  a  cutoff  frequency  which  is  dependent  upon  the 
microphone  diameter  and  the  Reynolds  number.  This  cutoff  frequency 
decreases  as  the  sensor  diameter  increases  because  the  spatial  averaging 
characteristics  or  larger  sensors  tends  to  cancel  out  much  of  the 
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convective  wavenumber  portion  of  the  wavenumber/ frequency  spectrum.  These 
examples  illustrate  that  the  overall  level  of  G  is  lowest  tor  the  larp-r 
measurement  sensors;  a  result  predicted  previously  bv  Girons  (oj.  Above 
the  cutoff  frequency,  the  spectral  level  decreases  approximately  is  * . 
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LLST  OF  FIGURES 


Schematic  description  of  boundary  layer  flow  over  a  flat  surface 
that  contains  acoustical  pressure  measurement  sensors.  The  wall 
is  assumed  to  be  rigid  and  u(-)/c  is  assumed  small,  where 
c  -  sound  speed. 

Non-Dimensionalized  bias  for  intensity  estimate  of  liquation  (17) 
as  a  function  of  kcAx  or  Ax/A. 

Wavenumber  response  function  for  circular  sensor.  The  experimental 
curve  is  from  Ref.  [14]. 

Schematic  description  of  the  wavenumber/ f requencv  spectrum,  where 

(a)  shows  the  ui  -  plane  and  (b)  shows  the  f>  -  kj  plane. 

Theoretical  wavenunoer/ f requency  spectra  based  on  the  models  of 
Maestrello  [17]  and  Chase  [18]. 

A  comparison  of  the  wall  pressure  spectrum  measured  by  Blake  [13] 
with  a  theoretical  prediction  based  on  liquation  ( A—  1 )  and  the 
Chase  [18]  model  for  i(kj .kj.w) . 

.  Example  calculations  of  G'jp(  „o)  for  various  values  of  the  r  low 
Reynolds  number  and  for  sensor  diameters:  (a)  13. 4o  nn, 

(b)  12.70  min,  (c)  0.35  mm,  and  (d)  3.18  ran. 
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FREQUENCY  (Hz) 

of  the  wal 1  pressure  spectrum  measured 
with  a  theoretical  prediction  based 
A— 1 )  and  the  Chase  [18]  model  tor 
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Figure  7.  Example  calculations  of  fl  {:.  )  for  various  values 
of  the  flow  Reynolds  number  and  for  sensor 
diameters:  (a)  _’).«(>  mm,  (li)  17.70  mm,  (e) 

h  .  19  n.n ,  and  (  d  )  5 .  !  S  i::m . 
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